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The correlation buildup and the formation dynamics of the shell structure in a spherically confined 
one-component plasma are studied. Using Langevin dynamics simulations the relaxation processes 
and characteristic time scales and their dependence on the pair interaction and dissipation in the 
plasma are investigated. While in systems with Coulomb interaction (e.g. trapped ions) in a 
harmonic confinement shell formation starts at the plasma edge and proceeds inward, this trend is 
significantly weakened for dusty plasmas with Yukawa interaction. With a suitable change of the 
confinement conditions the crystallization scenario can be externally controlled. 
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Crystallization of charged particles, predicted by 
Wigner seven decades ago, continues to stimulate re- 
search in many fields due to its relevance for astro- 
physics (e.g. white dwarf stars), basic many-body physics 
and potential applications in quantum computing, e.g. 
yy. Experimental realizations include electrons on a he- 
lium surface ion in traps [H electrons in quan- 
tum dots [4], dusty plasmas [Sj |6] and ultracold neu- 
tral plasmas [71 [8]. For crystallization to be possible 
in thermodynamic equilibrium in a macroscopic three- 
dimensional one-component plasma (OCP), the coupling 
parameter, T = q^/{av,sk^T) has to exceed a value of 
about 174, where Ows denotes the Wigner-Seitz radius, q 
the charge, T the temperature and Boltzmann's con- 
stant, whereas additional conditions have to be fulfilled 
in two-component plasmas (TCP) 9 . Recently, crys- 
tallization of spherically trapped dust particles has been 
achieved [lOj which revealed close similarities to ion crys- 
tals in traps [U [3] with the main difference being the 
screening of the Coulomb interaction in the former case 
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While the stuctural properties of the crystals are well 
understood, e.g. [HI US], much less is known on the dy- 
namics of their formation. Murillo showed [TJ that a 
neutral TCP produced by rapid ionization of an atomic 
gas will not crystallize because the correlation buildup is 
accompanied by heating |15| which limits F. Then Pohl 
et al. demonstrated that crystallization can be achieved 
if the expanding plasma is laser-cooled [T5], which still 
has not been realized experimentally [71 [Sj. An inter- 
esting prediction of [16] was that spherical crystal shells 
start to form in the cluster core. It is an open ques- 
tion whether this is a general crystallization scenario in 
trapped plasmas since, so far, no investigations on the 
crystal formation dynamics in spherically trapped ions 
and dusty plasmas have been performed. 

The goal of this Letter is, therefore, to perform such 
an analysis for spherically confined dusty plasmas. We 
study in detail the time-dependence of crystallization by 



simulating an experimental cooling process from a weakly 
correlated finite dust cloud towards strong coupling. We 
show that the formation of spatial correlations proceeds 
in a sequence of stages and present results for the char- 
acteristic time scales. Further, the dependence of the 
dynamics on screening and dissipation is explored. We 
predict that the onset of shell formation is typically at 
the cluster edge but the order of appearence of the in- 
ner shells can be controlled by suitable variation of the 
confinement. Finally, when the core region is made un- 
accessible to the plasma crystallization can be initiated 
in the center. 

Model and simulation idea. We consider N identical 
particles with mass m and charge q interacting through 
a Yukawa pair potential [Hj (f){r) — q^e~'"'/r in an ex- 
ternal confinement V{r). The effective range of (f){r) is 
determined by the screening parameter k. The dynamics 
of our system is described by the Hamiltonian 

i—l i—1 i^'j 

U(vi,....v„) 

With a harmonic confinement V{r) = mujQr^ /2 this 
model accurately describes the properties of spherical 
dust balls observed in recent experiments |llj and is 
equally applicable to spherically trapped ions in the limit 
K = 0. Below we use dimensionless units with the char- 
acterstic length and energy scales a = (^^/toWq)^/'^ and 
= q^/a. 

The ambient neutral gas in dusty plasma experiments 
is accounted for by an additional damping term and a 
fluctuating force in the (Langevin) equation of motion of 
the i-th particle 

mvi = -ViUiyi, . . . ,rjv) - mvi, + i^{t). (2) 

The friction coefficient v and the Gaussian noise 
^i{t) are related by the fluctuation-dissipation theorem 
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^ff (t)fj'(t')) = 2miykBT*S,jSapS{t - t'), where a,/3 G 
{x, y, z} and i, j G {1, . . . , N}. In Ref. [TB] this method 
was used to describe the cooling effect of the laser. 

We consider the following scenario to study the buildup 
of correlations: We start from a weakly correlated steady 
state of N trapped dust particles characterized by Fi = 

/{ak-QT{) = 0.2 which can be prepared e.g. by contin- 
uous laser-heating, which has been successfully applied 
in experiments on 2D layers [T7]. At the time t — Q the 
laser is turned off and the particles begin to slow down at 
a rate determined by the friction coefficient v, approach- 
ing a new equilibrium at the neutral gas temperature Tn 
corresponding to strong coupling with V-^ = 125. This 
scenario allows us to study the correlation buildup in 
a well defined manner without introducing spatial inho- 
mogeneities. The chosen Tn is of the order of recently 
measured normal mode temperatures [18j . Due to the 
external trap the system is inhomogeneous but isotropic, 
i.e. the mean density will only depend on the distance r 
from the trap center. To reduce the effect of fluctuations 
we perform several hundred to 1200 runs with random 
initial conditions over which we average. 
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FIG. 1: (Color online) Time dependence of the kinetic, po- 
tential and interaction energy for A'^ = 400. For the sake 
of clarity the graphs for different parameters are shifted by 
A(a;ot) = 5 (left) and l^(yt) = 1.5 (right). Left: Influence of 
screening for v jut^ = 0.2 and na = 0.6, 0.4, 0.2, (from left 
to right). Right: Influence of the damping rate for Ka — 0.6 
and v/u>o = 0.2, 0.5, 1, 2.5. Note the different scaling of the 
time axes. Potential and interaction energy are normahzed to 
the equihbrium energies at Fi = 0.2. 

Cooling towards strong coupling. Consider first the 
evolution of the coupling parameter T(t) which we 
compute from an instantaneous temperature k^Tit) — 
2E],init)/3N [19]. Fig. [l] shows that F(t) increases con- 
tinuously reaching the value 125 within uot ^ 35 |20j . 
The increase of F(t) is accompanied by weak oscillations 
and is only marginally affected by k, cf. left column. In 
contrast, the influence of friction is more apparent: an 
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FIG. 2: (Color online) Left: Instantaneous oscillation fre- 
quency of Epot (cf. Fig. [l| for various k at u/uio = 0.1. 
Horizontal lines denote the mean frequency in the interval 
35 < ujot < 50. Right: Mean frequency (corrected for friction, 
see text) compared to the analytical expression of Ref. |21[ . 



increase of v leads to an increase of the modulation am- 
plitude, cf. right column of Fig. [l] For strong friction 
the initial growth follows T{t) cx exp(2z^t), which is the 
expected behavior for a free particle subject to friction 
(baUistic regime). 

Fig. [l] also shows a non-trivial dynamics of the con- 
finement (Epot) and interaction energy (E-mt) during the 
crystallization. While Ep^t decreases, due to compres- 
sion of the cluster, Ei^t inreases due to the formation 
of correlations. The relative gain (loss) of interaction 
(confinement) energy increases with screening. E'pot and 
i?int exhibit much more pronounced oscillatory modu- 
lations than r(<) which are determined by dissipation. 
While for ly/ujQ = 0.2 only small oscillations occur, for 
v/uq = 0.5 an overshooting of i?pot is observed which 
reaches its maximum for v/too — 1. Upon further in- 
crease of I' the oscillations vanish. 

The origin of the oscillations is easy to understand. 
When the heating is turned off the amplitude of the ran- 
dom force is reduced by A/o = y/viTi — y/vTn, giving 
rise to a rapid radial contraction of the cluster which ex- 
cites a monopole oscillation. An increase of v leads to 
a faster loss of kinetic energy and a stronger contrac- 
tion, explaining the larger oscillation amplitude. For 
v ^ 2wo the oscillation is overdamped and the ampli- 
tude decreases, whereas for v <ti luq, A/q is small and 
the system smoothly evolves from one equilibrium state 
to another. Therefore, there exists a maximum in the 
oscillation amplitude observed at v « loq. Let us now 
analyze the oscillation frequency. To this end we com- 
pute an instantaneous frequency uj*([ti + <i+i]/2) from 
two successive minima or maxima of the potential en- 
ergy at ti and ti+i, see left part of Fig. [2] After a few 
cycles the frequency saturates and the oscillations corre- 
spond to a damped normal mode of the new equilibrium 
state. The intrinsic normal mode frequency uj (of the 
dissipationless system), right part of Fig. [2] is computed 
via UJ* = {uj'^ — i/^/4)^/^ [21]. For Coulomb interaction 
UJ agrees with the breathing frequency ujhr = V^ujq |12j . 
whereas in the case of finite screening it depends on k. 
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In the right part of Fig. [2] we also display an analytical 
result for Wbr(K^o) derived for a homogeneous Yukawa 
sphere |21j which is accurate at low screening. The nor- 
malized radius kRq is computed from a mean-field the- 
ory [21 m], see Eq. ^ below. 

Time- dependent density profile. The evolution of the 
radial density profile is shown in Fig. [3] In the initial 
weakly coupled state (Fi = 0.2) the density is monoton- 
ically decaying and is well described by the Boltzmann 
factor, noff(r) oc exp[—Vcs{r)/kBTi], where Vcs denotes 
the sum of confinement and mean field potential [25] . 
As the kinetic energy drops the system approaches the 
strong coupling regime with its characteristic shell struc- 
ture. In the case of Coulomb interaction (a, d) the first 
shell appears around the time ujQt w 10, at the cluster 
boundary. Only upon further cooling [increase of F(i)] 
shells form one after another in the direction of the trap 
center after almost constant time intervals, cf. dashed 
line in Fig. |3]d. The situation is different in the case of 
Yukawa interaction. While shell formation again starts 
at the edge, inner shells form more rapidly, cf. dashed 
line in Fig.[3je, and at kg = 2, (Fig.[3jf), the inner shells 
emerge almost simultaneously. 

This sequence of shell formation is in striking contrast 
to the one observed in expanding laser-cooled plasmas 
[IB] where shells emerge in the center, which was at- 
tributed to an increased density in the core. In fact, 
the mean density profile n(r) allows one to define a local 
coupling parameter T{r) ~ fi(r)^/'^, and it is tempting to 
exepect shell formation to start at a radius where T{r) 
has its maximum. Our results allow us to verify this 
hypothesis. During the initial phase, < LOot ^ 10, h 
evolves from the Boltzmann factor, ncff(r), to a (still) 
monotonous profile, cf. Fig. |3]a-c, which is very well 
described by the zero-temperature mean field result no 
given by [23] 

no{r, N) = [AV{r) - KV(r) + /W, (3) 

where fi{N) is the chemical potential, except for a slightly 
smoother decay at the edge. Note that Eq. (|3| exhibits a 
finite density step An(i?o) at some maximum radius Rq 
[noir, N) = for r > Rq] which emerges in our simula- 
tions rapidly, within toot ^ 10. 

In the Coulomb case [Fig.[3ja], n{r) is almost constant 
for LUot < 10, with a slight decay towards the edge, in 
agreement with Eq. ([3| which predicts a constant den- 
sity for r < Rq = 10.6. Nevertheless, shells appear at 
very different moments starting at the edge, where the 
density is smallest. For kq — 2, the mean density de- 
creases even stronger, f (i?)/f (0) w (1/2)^/3^ Fig.|3]c, 
but even here shells form at the edge first, in contrast to 
the above expectation. In fact, just prior to formation 
of the first shell around toot ~ 10 the coupling parameter 
r(t) approaches 10, cf. Fig. [l| where the mean-field de- 
scription fails and correlations become important. With 
increasing T(t) the discrete nature of the particles begins 
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FIG. 3: (Color online) Right: Evolution of density profile for 
TV = 1200 with v/luo ~ 0.2, black dashed line connects the 
approximate times and positions at which the shells emerge. 
Left: Snapshots of n(r) at ujot = 0, 15, 30 (oscillations grow 
with time) together with solution of Eq. ([3|, "MF", black 
dots. From top to bottom row: na — 0; 0.6; 2 



to manifest itself leading to formation of a (correlation) 
"hole" around each particle which cannot be occupied by 
others. This separation of particles in radial (and tan- 
gential) direction and an overall expansion of the cluster 
cause an increase of i?pot proportional to the number of 
particles in the edge layer, A'odgo ~ An{R). The system 
reduces this energy by spontaneously restructuring such 
that particles from the edge are accumulated at a smaller 
distance - the outer shell forms. The formation of the in- 
ner shells is triggered by the continuing increase of F in 
the center. The substantial acceleration of shell forma- 
tion with increasing k is explained by the inward force 
exerted by the particles on the outer shell [23] which is 
neghgible in the Coulomb case (Faraday cage effect). 

Effect of the confinement potential. The different crys- 
tallization behavior observed in the expanding neutral 
plasma can be traced to a different confinement potential. 
There, the ions are confined by the mean-field potential 
of the electrons [16] which is Coulomb-like (except for the 
center), i.e. V{r) ^ r~^. Consequently the mean density 
profile, prior to shell formation, exhibits a drastic in- 
crease towards the center, fi{r) ^ r^^, cf. Eq. and a 
vanishing density step An(i?) 0. The latter arises from 
vanishing of V{r) for r — > oo, i.e. ions with a finite kintic 
energy cannot be confined; instead of accumulating par- 
ticles in a shell, the system expands. The strong increase 
of n(r) towards the center then explains the observed 
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FIG. 4: (Color online) Equilibrium density profiles for dif- 
ferent trap potentials and various temperatures: a) quartic 
confinement V4, b) linear confinement Vi. c) and d): har- 
monic potenital with blocked core, V2, r^/a — 2. Units are, 
ri = q/y/ci, TA = \f(FJcA and ri,4 = g^/(?-i,4 fca). 



dynamics. We now further verify the governing role of 
the confinement potential for the dynamics of shell for- 
mation. To this end we analyze the equilibrium density 
profile |25j for different temperatures by performing ther- 
modynamic Monte Carlo simulations for a finite number 
of particles trapped by different potentials. According 
to Eq. ([3]) a potential Vir) cx r" yields a mean density 
n(r) cx r"^^, for k = 0. Fig.|4^. shows results for a quar- 
tic potential, V4(r) = C4r'*/4, where n^ir^ oc . Here a 
very large density step is formed which, together with the 
radial density increase, further enhances the shell forma- 
tion from the edge, as for the harmonic confinement, cf. 
Fig. |3^.,d. Next, consider a (regularized) linear confine- 
ment, Vi(r) = ci[rj +r^]-'^/^, for which n(r) cx r^^ . While 
again the first density maximum emerges at the bound- 
ary, strong modulations of n(r) appear near r = at 
low temperatures. Fig. |4]3, i.e. crystallization proceeds 
simultaneously from the edge and from the center. 

Finally, we further modify the confinement by making 
the central part of the trap, < r < r^, unaccessible for 
the particles by using an infinite wall at r = together 
with a shifted harmonic potential V2(r) (r — r^)^, cf. 
inset of Fig. |4}i. The results for the cases of Coulomb 
and Yukawa interaction (ko = 1) are strikingly different, 
cf. Fig. |4j:,d. While in the former case again shell forma- 
tion starts at the edge, in the latter the first shell clearly 
emerges in the core. The reason is that for Coulomb 
interaction, particles at r = r^, experience almost no 
radial force, whereas for Yukawa interaction charges lo- 
cated outside do produce an inward force [25] causing 
strong particle accumulation at the wall. Not only does 



this allow to reverse the spatial crystallization dynamics, 
this is also a situation where qualitatively different be- 
havior should be observable for spherically trapped ions 
(Coulomb interaction) and dusty plasmas (screened in- 
teraction) . 

In Summary, we have studied the transition from a 
weakly coupled to a strongly coupled state in a spheri- 
cally trapped dusty plasma in a scenario which can be 
realized experimentally. The initial relaxation phase, 
LJat ^ 10, is characterized by formation of a density 
step An at the edge and, for low friction, excitation of 
a breathing oscillation with a k— dependent frequency. 
Around the time co^t ~ 10 shell formation starts at the 
plasma edge which is a correlation effect arising from the 
finite density step. For typical dusty plasma experiments 
with Kfl « 0.6 [TT] and harmonic confinement, inner shells 
are formed one by one within almost constant time inter- 
vals ol wqI ~ 4. Furthermore, the crystallization dynam- 
ics can, to a large extent, be controlled by the shape of 
the confinement potential. In particular, by blocking the 
central part of the trap crystallization can be initiated in 
the center. 
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